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Abstract. Satellite CMB anisotropy missions, such as WMAP and Planck, and 
also the new generation of balloon-borne and ground experiments, make use of 
complex multi-frequency instruments at the focus of a meter class telescope to 
allow the joint study of CMB and foreground anisotropies, necessary for a high 
quality component separation. Between ~ 70 GHz and ~ 300 GHz, where fore- 
ground contamination is minimum, it is extremely important to reach the best 
trade-off between the improvement of the angular resolution, necessary to mea- 
sure the high order acoustic peaks of CMB anisotropy, and the minimization of 
the straylight contamination mainly due to the Galactic emission. This is one of 
the most critical systematic effects at large and intermediate angular scales (i.e. 
at multipoles £ less than ~ 100) and constists in unwanted radiation entering the 
beam at large angles from the direction of the antenna boresight direction. We 
focus here, as a working case, on the 30 and 100 GHz channels of the Planck Low 
Frequency Instrument (LFI). By assuming the nominal Planck scanning strat- 
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egy, we evaluate the GSC introduced by the most relevant Galactic foreground 
components for a reference set of optical configurations, accurately simulated as 
described in Paper I. We show that it is possible to improve the angular resolu- 
tion of 5 — 7 % by keeping the overall peak-to-peak GSC below the level of few 
jj,K (and about 10 times smaller in terms of RMS). A comparison between the 
level of straylight introduced by the different Galactic components for different 
beam regions (intermediate and far sidelobes) is presented. Simple approximate 
relations giving the RMS and peak-to-peak levels of the GSC for the intermedi- 
ate pattern and far sidelobes as functions of the corresponding contributions to 
the integrated antenna pattern are provided. We compare the results obtained at 
100 GHz with those at 30 GHz, where GSC is more critical. Finally, for some 
reference cases we compare the results based on Galactic foreground templates 
derived from radio and IR surveys with those based on WMAP maps including 
CMB and extragalactic source fluctuations. 

Key words. Cosmology: cosmic microwave background - Galaxy: general - Space 
vehicles - Telescopes - Methods: data analysis. 

1. Introduction 

After the detection of cosmic microwave background (CMB) anisotropies at few degree 
scales by COBE/DMR (Smoot ct al. 1992, Bennett et al. 1996, Gorski ct al. 1996) 
and the recent balloon-borne and ground experiments (sec Bcrsanclli et al. 2002 and 
references therein for a review on the pre- WMAP observational status) at high sensitivity 
and resolution on limited sky regions probing a universe model with f2 tot ~ 1 (see e.g. 
Netterfield et al. 2002, Stompor et al. 2001, Pryke et al. 2002, and references therein) 
the NASA space mission WMAP (Wilkinson Microwave Anisotropy Probe, sec Bennett 
et al. 2003a) derived the CMB anisotropy angular power spectrum with unprecedent 
sensitivity and reliability (Hinshaw et al. 2003b, Kogut et al. 2003) and improved the 
accuracy in the determination of the most important cosmological parameters (Spergel 
et al. 2003). 

Future fundamental progresses in CMB anisotropy and polarization will be based on 
the Planck mission by ESA 1 (Bersanelli et al. 1996, Tauber 2000, Villa et al. 2003), 
planned to be launched in the year 2007. 

In particular, the Low Frequency Instrument (LFI, Mandolesi et al. 1998; see also 
Mandolesi et al. 2002) and the High Frequency Instrument (HFI, Puget et al. 1998; 
see also Lamarre et al. 2002) on-board Planck will cover together a wide frequency 
range (30-900 GHz) which should significantly improve the accuracy of the subtraction 
of foreground contamination from the primordial CMB anisotropy, providing at the same 

Send offprint requests to: C. Burigana 
1 http://astro.estec.esa.nl/Planck/ 
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time a gold mine of cosmological as well as astrophysical information (see e.g. De Zotti 
et al. 1999 and references therein). 

To fully reach these scientific goals, great attention has to be devoted to properly 
reduce and/or subtract all the possible systematic effects. 

The effect of optical distortions, both in the main beam and in the near and far 
sidclobes, has been widely recognized as one of the most critical systematics both in 
balloon experiments and in space missions (Page et al. 2003, Barnes et al. 2003). 

The systematic effect introduced by main beam distortions (Burigana et al. 1998) can 
be in part reduced by adopting aplanatic configurations for the primary mirror (Villa et 
al. 1998, Mandolesi et al. 2000a, Villa et al. 2002) and its effect on the CMB anisotropy 
power spectrum recovery, relevant at the multipoles of the acoustic peaks, can be partially 
removed through dedicated deconvolution codes (Arnau et al. 2002) provided that the 
main beam shape can be accurately reconstructed in flight (Burigana et al. 2002). 

The requirement on the rejection of unwanted radiation coming from directions far 
from the optical axis (straylight) is stringent for Planck and does not pertain only the 
telescope itself, but the entire optical system, including solar panels, shielding, thermal 
stability and focal assembly components. The variations of the spurious straylight signal 
introduce contaminations in the anisotropy measurements. The removal of this effect in 
data analysis is in principle much more complicated than the subtraction of main beam 
distortion effect. This is due to the difficulty to accurately known the "real" antenna 
pattern (with both at ground and in flight reconstruction) at very low response levels. 

The antenna response features far from the beam centre (sidclobes) are determined 
largely by diffraction and scattering from the edges of the mirrors and from nearby 
supporting structures. Therefore, they can be reduced by decreasing the illumination at 
the edge of the primary, i.e. increasing the edge taper (ET if expressed in dB; the linear 
edge taper, LET, is 10~ ET / 10 ), defined as the ratio of the power per unit area incident 
on the centre of the mirror to that incident on the edge. Of course, the higher is edge 
taper, the lower is the sidelobe level and the straylight contamination. On the other hand, 
increasing the edge taper has a negative impact on the angular resolution for a fixed size 
of the primary mirror (see e.g. Mandolesi et al. 2000b). 

In the "cosmological window" between ~ 70 GHz and <~ 300 GHz, where foreground 
contamination is minimum, it is extremely important to reach the best trade-off between 
the improvement of the angular resolution, necessary to measure the high order acoustic 
peaks of CMB anisotropy, and the minimization of the straylight contamination due to 
the Galactic emission (GSC, Galaxy Straylight Contamination), one of the most crit- 
ical systematic effects, most relevant at large and intermediate angular scales (i.e. at 
multipoles £ less than « 100). 

In this work we focus, as a working case, on the 100 GHz channels of Planck Low 
Frequency Instrument, although the methods and the basic results described here can 
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be applied also to other Planck frequency channels and to different CMB anisotropy 
experiments. We shall also compare our results with simulated data at 30 GHz. 

We will use here the detailed analysis on optical computations and the main optical 
results concerning the full antenna pattern response given in Sandri et al. 2003 (here- 
after Paper I). We present extensive simulations of the GSC due to Galactic foreground 
components relevant at 100 GHz and, at this purpose, we have made use of a wide set of 
simulated optical configurations in order to find the best compromise between resolution 
and GSC therefore defining the maximum range of multipoles accessible to the considered 
frequency channel and the level of straylight signal affecting the data. 

In Sect. 2 we briefly describe the basic recipes to simulate Planck observations 
(see also Appendix A), the adopted optical input, and maps of Galactic components. 
The results of the simulated straylight contamination are described in Sect. 3, while the 
comparison between the results obtained for different antenna patterns and foreground 
components is presented in Sect. 4 (and in Appendix B for the comparison with a stray- 
light simulation at 30 GHz). Finally, we discuss the results and draw our main conclusions 
in Sect. 5. 

2. Simulations 

The selected orbit for Planck is a Lissajous orbit around the Lagrangian point L2 of 
the Sun-Earth system (see e.g. Mandolesi et al. 1998). The spacecraft spins at 1 r.p.m. 
and the field of view of the two instruments (LFI/HFI) is about 10° x 10° centered at 
the telescope optical axis (the so-called telescope line of sight, LOS) at a given angle 
a from the spin-axis direction, given by a unit vector, s, chosen to be pointed in the 
opposite direction with respect to the Sun. In this work we consider values of a ~ 85°, 
as adopted for the baseline scanning strategy. The spin axis will be kept parallel to the 
Sun-spacecraft direction and repointed by ~ 2.5' every ~ 1 hour (baseline scanning 
strategy). Hence Planck will trace large circles in the sky. A precession of the spin-axis 
with a period, P, of ~ 6 months at a given angle (3 ~ 10° about an axis, /, parallel to 
the Sun-spacecraft direction (and outward the Sun) and shifted of ~ 2.5' every ~ 1 hour, 
may be included in the scanning strategy, possibly with a modulation of the speed of the 
precession in order to optimize data transmission (Bernard et al. 2002). Although the 
scanning strategy could be changed, the GSC pattern, peak-to-peak, and angular power 
spectrum are very weakly dependent on the details of these proposed scanning strategies 
(Burigana et al. 2000). 

The code we have implemented for simulating Planck observations for a wide set of 
scanning strategies is described in detail in Burigana et al. (1997, 1998) and in Maino 
et al. (1999). In this study we exploit the baseline scanning strategy and simply assume 
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Planck located in L2. Wc do not consider the Planck Lissajous orbit around L2 
because its effects are negligible in this context. 

We compute the convolutions between the antenna pattern response and the sky 
signal as described in Burigana et al. (2001) by working at ~ 1° or ~ 7' resolution when 
considering the far or the intermediate sidelobes and by considering spin-axis shifts of 
~ 2° every two days and 180 samplings per scan circle. In fact, the effects of pattern 
features we want to study here occur at ~ degree or larger scales and, also, the wide set 
of optical simulations for the full antenna pattern is available at ~ degree resolution (see 
Paper I for details on computation time of optical simulations). 

With respect to the reference frames described in Burigana et al. (2001), following 
the recent developments in optimizing the polarization properties of LFI main beams 
(see Paper I), the conversion between the standard Cartesian telescope frame x,y, z and 
the beam frame Xbf,ybf,Zbf requires a further angle ipB other than the standard polar 
coordinates 9b and <j>B defining the colatitude and the longitude of the main beam centre 
direction in the telescope frame. Appendix A provides the transformation rules between 
the telescope frame and the beam frame, as well as the definition of the reference frames 
adopted in this work. 

The orientation of these frames as the satellite moves is implemented in the code. For 
each integration time, we determine the orientations in the sky of the telescope frame 
and of the beam frame, thus performing a direct convolution between the full pattern 
response and the sky signal for the desired number of maps, simultaneously. 

2.1. Optical inputs 

A detailed discussion of the optical simulation method and results is presented in Paper I 
(see Tables 2 and 3 of Paper I for the main properties of the adopted antenna patterns) . 
We briefly summarize here the most relevant aspect in this context. 

Several full beam patterns have been simulated for different designs of two Planck 
LFI feedhorns at 100 GHz differently located on the focal plane unit: three models (9A, 
9B, 9C) for LFI9 and four models (4A, 4B, 4C, 4D) for LFI4. Different values of ET will 
reflect in significant differences in the level of the GSC in the TOD, while different feed 
designs with the same ET will produce small, but not negligible, differences. 

The details of the antenna pattern response, computed as described in Paper I as 
functions of the two standard polar coordinates 6b f, 4>bf in the beam frame, depend also 
on the optical contributions considered in the analysis. For some representative cases, 
we will compare the results of our simulations of GSC by adopting optical computations 
taking into account the first and second or the first, second and third order optical 
interactions (see Paper I). 
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In the beam frame we identify two different angular regions relevant for the straylight 
analysis: the intermediate pattern and the far sidelobes, respectively defined as the region 
between 1.2° and 5° from the beam centre direction and at angles large than 5° from the 
beam centre direction. The definition of these "reference" angles is, of course, somewhat 
arbitrary. They roughly separate angular regions where significant pattern variations 
occur on sub-degree scales from those where they occur on scales of few degrees or larger. 
The values adopted here allow a direct comparison with a previous analysis (Burigana 
et al. 2001). 



2.2. Maps of Galactic components 

In the cosmological window (70-300 GHz) the CMB is clearly the dominant component; 
on the other hand the Galactic emission is still relevant at low and middle Galactic 
latitudes. The templates adopted here are similar to those described in Maino ct al. 
(2002) and Paladini et al. (2003). 

At 100 GHz the most relevant Galactic component is the thermal dust emission. We 
adopted here a template obtained by extrapolating the maps by Schlegel et al. (1998) 
which combines IRAS and DIRBE data, assuming a grey-body spectrum (expressed in 
antenna temperature), 

v f3+1 hv 

T A ,dust{v) oc , v = — (1) 

C t K-t-dust 

with uniform temperature T^ ust = 18 K and emissivity (3 = 2. 

In order to simulate the free-free contribution we assume, somewhat arbitrarily, that 
it is perfectly correlated with the dust itself, i.e. that it has the same spatial distribution. 
Its antenna temperature scales with frequency as Taj/ oc i/~@ff , with [3ff — —2.1. The 
relative amplitude of dust and free- free emission is assumed to be a factor of 3 at 100 GHz 
(De Zotti et al. 1999). Thus, we produce a single map of "thermal" emission from dust 
plus diffuse free-free emission (see Fig. 1) with a spectrum described by: 

\@ff , T A4ust (v) 
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The synchrotron emission template is the 408 MHz map of Haslam et al. (1982), 
available at a resolution of 0.85°, extrapolated to the considered frequencies (sec Fig. 2) 
assuming a uniform spectral index /3 syn = —2.9 in antenna temperature. No attempt is 
made here to add small scale fluctuations, since the effects of pattern features on which 
this work is focused occur at ~ degree or larger scales; for the same reason, the fact the 
original template includes a convolution with a beam with the given resolution is not a 
concern. 

Localized free-free emission dominates the signal on wide areas in the Galactic plane 
at least over the frequency range from 30 to 100 GHz. The Synthetic Catalog at 2.7 GHz 
produced by Paladini et al. (2003) provides a rich information on compact Galactic HII 
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dust + diffuse free-free free-free = 1/3 dust 
log T ant (mK) 




-3.3 — ■ — 1.2 



Fig. 1. Map of the Galactic dust emission plus diffuse free- free emission adopted in this 
study (see also the text). 



diffuse synchrotron 
log T ant (mK) 




-3.58 



-1.11 



Fig. 2. Map of the Galactic diffuse synchrotron emission adopted in this study (see also 
the text). 



regions that has been used in this work to generate a map of free- free emission from these 
sources (see Fig. 3). A spectral index a = —0.1 in flux (—2.1 in antenna temperature), as 
in the case of thermal bremsstrahlung emission in a thin plasma, has been adopted here to 
extrapolate the signal from 2.7 GHz to the considered Planck frequency channels. With 
respect to the simulated map reported by Paladini et al. (2003), we have implemented 
here a code which simulates the contribution of each source in the Synthetic Catalog 
to each map pixel without applying at this step the convolution with the beam, since 
the convolution with the intermediate and far antenna pattern has been subsequently 
applied as described in the first part of this section. 
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HII regions 
log T ant (mK) 




Fig. 3. Map of the free-free emission from compact Galactic HII regions adopted in this 
study (see also the text). 

All maps have been projected into the HEALPix scheme 2 (Hierarchical Equal Area 
and IsoLatitude Pixelization of the Sphere) by Gorski et al. (1999). 

While this work was nearly completed, the 1 year data products 3 from the WMAP 
satellite have become available. In a few representative cases, we have repeated the stray- 
light analysis by adopting the WMAP frequency maps at 33 and 94 GHz, including in 
the straylight evaluation also the minor contributions to straylight signal from CMB and 
extragalactic source fluctuations 4 . Note that WMAP maps are, of course, convolved 
with the corresponding beam patterns and include the effect of main beam distortions 
and straylight contamination as well as the instrumental noise and the effect of other 
sistematics not subtracted in the data analysis; on the other hand, these effects are sig- 
nificantly smaller than the signal (Hinshaw et al. 2003a) and, since this analysis is mainly 
contributed by signal variations on degree or larger angular scales, they can be neglected. 

3. Simulation results 

The main output of our simulation code are the time ordered data (TOD) of the signals 
entering the intermediate and far pattern. 

In the TOD of each scan circle, two prominent maxima typically appear. These are 
related, for the intermediate pattern, to the two crossings of the Galactic plane of the 
telescope field of view and to the crossings of the Galactic plane of the main spillover 

2 http:/ /www.eso.org/science/healpix/ 

3 http://lambda.gsfc.nasa.gov 

4 Note that these WMAP maps do not include the monopole, while our Galactic component 
maps include the corresponding monopoles. The comparison between the average straylight 
signals will be therefore only indicative. Clearly, anisotropy experiments are not sensitive to the 
monopole, directly subtracted in the data. 
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in the case of the far pattern. As already recognized by Burigana et al. (2001), these 
maxima are only slighlty shifted with respect to the maxima of the signal entering the 
main beam in the case of the intermediate pattern and shifted of about 90° in the case 
of the far pattern, as direct consequence of the pattern shape. 

The typical signal level is determined by the sky signal and the fraction, /% = 
100 J n JdQ/ J 47r JdQ, of integrated antenna response in the considered antenna pattern 
region, £t, reported in Table 3 of Paper I. The ratio between the fraction of the in- 
tegrated antenna response in the far pattern and the one in the intermediate pattern 
provides only a rough upper limit to the ratio of straylight pcak-to-peak signal in these 
two pattern regions. In fact, while an extended very bright Galactic region could quite 
easily fill the (relatively small) solid angle subtended by the intermediate pattern, the 
whole (quite large) solid angle subtcnted by the main spillover and by the other relevant 
far pattern features can not be easily filled by signals coming all simultaneously from 
very bright Galactic regions. Similarly, the fractional difference between the integrated 
antenna response in the far pattern computed by including or not the third order optical 
interaction provides an upper limit (less than about 10%) to the fractional underesti- 
mation of the GSC when the third order optical interaction is neglected. More accurate 
estimates require numerical simulations. 

In Figs. 4 and 5 we report the TOD corresponding to the straylight signal respectively 
from the far and intermediate pattern for the three Galactic components described in 
Sect. 2.2 in the case of the beam LFI9 9B computed by including the first, second and 
third order optical interactions. 

The TOD corresponding to the difference between the straylight signal obtained by 
including or neglecting the third order optical interaction are reported in Figs. 6 and 
7. Clearly, the difference is within ~ 5 % of the straylight signal, i.e. smaller than the 
above upper limit derived on the basis of simple optical considerations by about a factor 
two. Therefore, considering only the first and second order optical interactions does not 
introduce relevant loss of information in the optimization analysis of the optical design . 

By comparing Figs. 1, 2, and 3 with Figs. 4 and 5, note how the different angular 
distribution of the three considered Galactic components reflects in the pattern of the 
straylight signal TOD, more or less "diffuse" according to the considered component. 

Similar results are found for all the considered beam patterns, but with a peak-to- 
peak and RMS straylight signal significantly dependent on the adopted pattern shape. 
Instead of reporting the full simulation results, a concise comparison between the results 
obtained for the whole set of optical configurations is reported in the next section. 

5 On the contrary, as obvious, the best possible knowledge of the antenna pattern, including 
all possible effects, from those due dust and molecular contamination on mirror surfaces to those 
related to the mirror roughness and temperature behaviour, should to be taken into account in 
the final data analysis. 
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pi^ Axis Pointing (deg: 



Fig. 4. Synthetic view of the data stream in terms of decimal logarithm of the antenna 
temperature in /iK from all scan circles for the different Galactic components. The ecliptic 
coordinates properly refer here to the latitude of the telescope LOS (for graphic purposes, 
in this plot the range between —85° and —255° refer to the second half of each scan circle) 
and to the its longitude shift with respect to its initial direction, or equivalently to the 
shift of the spin axis pointing direction. We report here the straylight signal in the far 
sidelobes computed by including the first, second and third order optical interactions for 
the beam LFI9 9B (see also the text). 
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Fig. 5. The same as in Fig. 4, but for the signal from the intermediate pattern. In this 
case, Galactic HII regions produce a quite localized straylight contamination, null far 
from these sources and arbitrarily set to 10 _20 /xK for graphic purposes (as in Fig. 4, we 
report here the decimal logarithm of the antenna temperature in /xK; see also the text). 



4. Comparison between different antenna patterns and foreground 
components 

The statistical moments of the straylight signal TOD and its peak-to-peak value are 
reported in Tables 1-7 for our whole set of optical configurations at 100 GHz [in the 
tables, the data referring to the global straylight effect from intermediate pattern plus far 
sidelobes (I + F) are derived including also the contribution from the third order optical 
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Fig. 6. The same as in Fig. 4, but referring to the difference between the stray light 
signals computed by including or not the third order optical interaction (the antenna 
temperature in fiK in linear scale is here reported; see also the text). 
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Fig. 7. The same as in Fig. 6, but for the signal from the intermediate pattern (see also 
the text). 



interactions, when available]. Table 8 summarizes the basic information contained in the 
above tables as functions of the edge taper. 

As evident from the tables, the contamination from the far sidelobes is much more 
relevant than that from the intermediate pattern when the RMS of the TOD is considered; 
by looking at the peak-to-peak of the straylight signal, the contamination from the far 
sidelobes is larger by a factor of few units or comparable to that from the intermediate 
pattern. The larger impact of the far sidelobes with respect to the intermediate pattern 
is particularly remarkable for the case of the diffuse Galactic components, while it is 
less evident in the case the map of free-free emission from Galactic HII regions. We hnd 
in fact a quite general behaviour: more diffuse the component and more relevant the 
straylight contamination from far sidelobes with respect to the intermediate pattern. 

On the contrary, the skewness and kurtosis indices are larger for the straylight con- 
tamination from the intermediate pattern because they are more sensitive to localized 
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Table 1. Statistical moments and peak-to-peak of the simulated TOD (antenna tem- 
perature expressed in fiK) from intermediate and far sidelobes of the simulated beam 
pattern LFI9 9A at 100 GHz for different Galactic components. The row labelled with 
"1123" ("112") refers to the signal in the intermediate beam computed taking into ac- 
count the first, second and third (first and second) order optical interactions. The row 
labelled with "F123" ("F12") refers to to the signal in the far sidelobes computed taking 
into account the first, second and third (first and second) order optical interactions (see 
also the text). 

LFI9 9A 

BEAM SKEWNESS KURTOSIS 

REGION AVERAGE VARIANCE RMS PEAK-TO-PEAK INDEX INDEX 



DUST + DIFFUSE FREE-FREE EMISSION 



112 


2.27x10" 


-3 


5.81xl0~ 5 


7.63x10" 


-3 


2.61x10" 


-l 


9.45x10+° 


1.33xl0+ 2 


F123 


6.71x10" 


-2 


4.77xl0~ 3 


6.90x10" 


-2 


5.88x10" 


-l 


1.88x10+° 


4.92x10+° 


F12 


6.40x10" 


-2 


4.49xl0~ 3 


6.70x10" 


-2 


5.68x10" 


-l 


1.89x10+° 


4.90x10+° 


I + F 


6.93x10" 


-2 


4.83xl0~ 3 


6.95x10" 


-2 


5.88x10" 


-l 


1.81x10+° 


4.59x10+° 



DIFFUSE SYNCHROTRON EMISSION 



112 


2.10x10" 


-4 


1.05x10" 


■7 


3.24x10" 


-4 


4.66x10" 


-3 


6.59x10+° 


5.64xl0 +1 


F123 


5.87x10" 


-3 


1.65x10" 


•5 


4.07x10" 


-3 


2.32x10" 


-2 


1.27x10+° 


1.11x10+° 


F12 


5.60x10" 


-3 


1.55x10" 


■5 


3.94x10" 


-3 


2.24x10" 


-2 


1.29x10+° 


1.13x10+° 


I + F 


6.08x10" 


■3 


1.67x10" 


-5 


4.09x10" 


-3 


2.31x10" 


-2 


1.22x10+° 


9.44xl0 _1 


HII REGIONS 


112 


7.70x10" 


-5 


7.45x10" 


■7 


8.63x10" 


-4 


8.05x10" 


-2 


4.14xl0 +1 


2.67xl0+ 3 


F123 


2.13x10" 


-3 


7.14x10" 


■6 


2.67x10" 


-3 


2.67x10" 


-2 


1.93x10+° 


5.87x10+° 


F12 


2.02x10" 


-3 


6.69x10" 


■6 


2.59x10" 


-3 


2.61x10" 


-2 


1.94x10+° 


5.82x10+° 


I + F 


2.20x10" 


-3 


7.99x10" 


-6 


2.83x10" 


-3 


8.67x10" 


-2 


2.97x10+° 


3.27xl0 +1 



SUM OF THE ABOVE COMPONENTS 



112 


2.55x10" 


-3 


7.02xl0~ 5 


8.38x10" 


-3 


2.94x10" 


-l 


9.51x10+° 


1.37xl0+ 2 


F123 


7.51x10" 


-2 


5.71xl0" 3 


7.56x10" 


-2 


6.38x10" 


-i 


1.84x10+° 


4.66x10+° 


F12 


7.16x10" 


-2 


5.38xl0" 3 


7.33x10" 


-2 


6.16x10" 


-i 


1.85x10+° 


4.64x10+° 


I + F 


7.76x10" 


-2 


5.79xl0^ 3 


7.61x10" 


-2 


6.38x10" 


-l 


1.78x10+° 


4.34x10+° 



features. As expected, these indices are larger for the free-free emission from Galactic 
HII regions than for the more diffuse components. 

For a couple of representative cases, LFI9 9B and LFI4 4A, we evaluated the stray- 
light signal on the basis of the WMAP map at 94 GHz for the intermediate pattern 
and far sidelobes, respectively (see Tables 2 and 4). For the signal in the intermediate 
pattern, we find a straylight contamination larger by a factor ~ 1.6 in terms of RMS 
(but essentially unchanged in terms of peak-to-peak) . It is interesting to note that in the 
case of the WMAP map, since the 94 GHz channel has a low level of Galactic foreground 
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Table 2. The same as in Table 1, but for the indicated simulated beam pattern. Note 
how the straylight in the intermediate pattern 1123 the GSC is just smaller when the 
third order optical interaction is included (compare 1123 with 112). This is due to a small 
decrease (~ 0.4 %) of the integrated antenna pattern response in that pattern region when 
the third order optical interaction is taken into account because of the combination in 
amplitude and phase of the various contributions, producing an averall antenna response 
- different from a simple sum of the powers of the various contributions - that typically 
increases with the number of considered contributions, as occurs here in the case of far 
sidclobcs, but not necessarily does. For the straylight in the intermediate pattern we 
report also the result based on the WMAP map at 94 GHz, including all components. 
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Table 3. The same as in Table 1, but for the indicated simulated beam pattern. 
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contamination, the contribution to straylight from CMB fluctuations clearly appears (see 
Fig. 8) far from the Galactic plane, where the straylight signal is negligible in the case 
of the simulation carried out by using the three Galactic templates (on the contrary, 
this effect is not evident at 30 GHz, see Fig. B.l in Appendix B). If we set to zero the 
values of the TOD samples with signal smaller than ~ 0.15/iK (the maximum value of 
straylight far from the Galactic plane in the right panel of Fig. 8) we obtain a RMS stray- 
light of ~ 0.024^iK, quite close to the RMS straylight value of ~ 0.028^K obtained in the 
case of the sum of the three adopted Galactic templates (see Table 2) . This means that 
main contribution to the factor ~ 1.6 of discrepancy in the RMS of the straylight signal 
found by adopting the three Galactic templates or the WMAP map at 94 GHz is mainly 
produced by the CMB anisotropy. 

In the case of the straylight in the far sidelobes we find a straylight contamination 
larger only by a factor ~ 1.3 in the case of the WMAP map than in the case of the sum of 
the three adopted Galactic templates, in terms of both RMS and peak-to-peak. The above 
CMB anisotropy effect is not evident. On the contrary, the regions at low straylight signal 
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Table 4. The same as in Table 1, but for the indicated simulated beam pattern. For 
the straylight in the far sidelobes we report also the result based on the WMAP map at 
94 GHz, including all components. 
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in the case of the simulation based on the WMAP map show a pattern that, although 
with different signal values, is quite similar to that found in the same regions by using the 
Galactic templates. This indicates that the main difference found in this case is due to the 
adopted model of Galactic templates, in agreement with the idea that CMB anisotropics 
(and also extragalactic bright sources and fluctuations, minimum at about 100-200 GHz) 
contained in the WMAP map, can produce only minor contaminations, compared to the 
Galactic signal, because of their smaller power at large angular scales which implies that 
positive and negative fluctuations partially compensate inside the relatively large antenna 
pattern solid angles relevant for this kind of straylight contamination. 

By exploiting our set of optical configurations and taking into account the correction 
factor found on the basis of the WMAP map, we find a linear approximation describing 
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Table 5. The same as in Table 1, but for the indicated simulated beam pattern. 
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Fig. 8. The same as in Fig. 5 but for the sum of the three Galactic templates and for 
the WMAP map at 94 GHz where the contribution from CMB anisotropy appears (the 
antenna temperature in [iK in linear scale is here reported; see also the text). 

quite well the dependence of the RMS and the peak-to-peak values of the GSC at 100 GHz 
on the (per cent) fractional contribution, f%, to the integrated antenna pattern response 
from the considered pattern region (see also Fig. 9): 



peak-to-peak ~ 5.3/iK x f% ~ 175^K x LET 



(3) 



rms ~ 0.62^K x f % ~ 21^iK x LET , 



(4) 
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Table 6. The same as in Table 1, but for the indicated simulated beam pattern. 
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for the far sidelobes, and 

peak-to-peak ~ 28/jK x f% ~ 98^K x LET (5) 

rms ~ 0.86^K x f% ~ 3^K x LET , (6) 

for the intermediate pattern 6 , where eqs. (l)-(5) of Paper I are taken into account in 
the last equalities of these equations (the LET is here referred to a taper at an angle 
of 24° from the feed axis). The numerical coefficients in eq. (6) can be multiplied by 
~ 1.6 to include also the stray light contribution in the intermediate pattern from CMB 
fluctuations. 

Clearly, in the cosmological window the GSC contamination is less than that in the 
lowest and highest Planck frequency channels and so its impact on CMB power spec- 
trum recovery (see Burigana et al. 2001 for an analysis of GSC impact on the power 
spectrum recovery and Fourier decomposition of scan circle data at the lowest Planck 
frequency channel). On the other hand, the ultimate goal of Planck, and in general of 
future CMB anisotropy experiments after WMAP, is not only the power spectrum recov- 
ery but also a detailed imaging of the last scattering surface and a detailed study of the 
6 In these fits we adopt, for uniformity, the numbers found by considering the first and second 
order optical interactions, available for the all cases. 



Rease give a shorter version with: \authorrunning and/or \titilerunning prior to \maketitle 
Table 7. The same as in Table 1, but for the indicated simulated beam pattern. 
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Table 8. RMS and peak-to-peak of the simulated TOD (antenna temperature expressed 
in fiK) from intermediate, far sidelobes and both of them of the all simulated beam 
patterns at 100 GHz for sum of the three considered templates of Galactic components as 
functions of the edge taper. We neglect here the contribution from the third order optical 
interactions since only first and second order optical interactions have been evaluated for 
the all optical designs. 
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whole statistical information, cosmological and astrophysical, contained in the frequency 
maps. For this reason, the Planck requirement on the maximum acceptable level of 



Please give a shorter version with: \authorrunning and/or \titilerunning prior to \maketitl© 



. - - ■ □ 



O . O O . 2 O . 4- O - 6 O . S 1 . O 

Fig. 9. Dependence of the RMS (solid line for the intermediate pattern and triangles for 
the far sidelobes) and the peak-to-peak values (dashed line for the intermediate pattern 
and squares for the far sidelobes) of the GSC at 100 GHz on the (per cent) fractional 
contribution to the integrated antenna pattern response from the considered pattern 
region and corresponding linear approximations (dotted lines). 

systematic effect contamination is of few ^K, in order to avoid spurious signals at a level 
comparable with the Planck sensitivity. 

Among the set of analyzed optical configurations we can identify a subsample which 
reaches a good trade-off between the angular resolution (see Table 2 of Paper I) and the 
suppression of the straylight contamination. 

In the case of LFI9 the configuration 9B shows a global peak-to-peak (rms) GSC 
less than about 2/zK (0.2£iK) simultaneoulsy reaching an angular resolution (FWHM) of 
10.02'. The configuration 9A shows a lower GSC, but with a worse angular resolution 
(10.56'), while the configuration 9C allows to reach a resolution of 9.54' but introducing 
a relatively large GSC, being its peak-to-peak (rms) signal of about 7/j,K (0.7/xK). The 
worse resolution for LFI4 with respect to LFI9 is due to the different location of the 
corresponding feed in the focal surface (see Mandolcsi ct al. 2000b). On the other hand, 
our analysis identifies a well defined optimal optical configuration (4C) that allows to 
reach the best angular resolution (12.08') with the minimum GSC, i.e. a peak-to-peak 
(rms) signal less than about 1.5/zK (0.15/xK). 

In Appendix B we report also the results obtained for a simulation at 30 GHz. We 
find some qualitative differences between the GSC at these two frequencies, related to 
their optical behaviours and to the different role of the Galactic foregrounds. 

An accurate computation of the GSC on Planck polarization data could be in prin- 
ciple carried out by using the formalism described in Challinor et al. (2000). On the 
other hand, since no microwave polarization surveys are currently available, a detailed 
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computation based on current simulated templates (see e.g. Giardino et al. 2002) may 
provide only indicative results. A first order analysis aimed to obtain a robust upper limit 
on the RMS GSC in the Planck LFI polarization data can be obtained with a simple 
argument. The signals, (i = 1,2), in the four radiometers associated to a couple of 
LFI feeds symmetrically located on the focal plane can be expressed as: 

2/u = T + 6Th + (Q + <yQ; 1 )cos(20n) + (U + SU^an^i) 

2/12 = T + 5T S 12 - (Q + <5Q* 2 )cos(20ii) - (U + <5C/ 1 s 2 )sin(20ii) 

27 21 = T + 6T& + (Q + <5Q^)cos(20 21 ) + (U + <5C/ 2 5 1 )sin(20 21 ) 

2/ 22 = T + <5T| 2 - (Q + 5Q s 22 )cos(2cf> 21 ) - (U + SU& )sin(20 21 ) . 

Here foi = </>n + 7r/4, 0n is the angle between the axis Xbf of the 6eam frame corre- 
sponding to the first feed and the direction of the parallel in the considered pointing 
direction, T, Q and U are the unpolarized signal and the Stokes parameters in the 
main beam, and the terms ST-j , 5Q*j , 5Q*j account for the spurious contributions due 
to the Galactic straylight. For each pointing direction, they are different for the four 
radiometers because of the different level and orientation of the corresponding antenna 
patterns. The measure of the Stokes parameters Q and U is obtained by combining 
the signals in the four radiometers: Q = (In — Ji2)cos(20n) — (I 21 — 7 22 )sin(2(/>ii), 
U = (I21 — /22)cos(20n) + (In — /22)sin(20n). The Galactic synchrotron emission is 
partially polarized (~ 30 %), while the Galactic free-free and thermal dust emission, 
most relevant at ^ > 50 GHz, are only weakly polarized (a few %). Therefore, at least 
in the cosmological window, a Galactic straylight contamination relevant for the polar- 
ization measure mainly derives from the differences between the temperature straylight 
signals, in each pairs of radiometers associated to the same feed. Assuming a typical 
value ~ 1/2 for cos(20n) and sin(20n), we find a RMS GSC on Q and U similar to that 
found above for T [see eqs. (4) and (6)], to be considered as a pessimistic upper limit, 
corresponding to a difference of a factor ~ 2 in the temperature straylight signals in each 
pair of radiometers. 

We then conclude that, at least in terms of RMS and in the cosmological window, 
keeping at very low level the GSC in the temperature data directly assures an adequate 
suppression of the GSC in the polarization data. 

5. Conclusions 

Satellite CMB anisotropy missions, such as WMAP and Planck, and also the new gen- 
eration of balloon-borne and ground experiments, make use of complex multi- frequency 
instruments at the focal surface of a meter class telescope to allow the joint study of CMB 
and foreground anisotropics, necessary for a high quality component separation. In the 
so-called "cosmological window", between ~ 70 GHz and ~ 300 GHz, where foreground 
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contamination is minimum, it is extremely important to reach the best angular resolu- 
tion (necessary to measure the high order acoustic peaks of CMB anisotropy) achievable 
keeping at the same time the straylight contamination at acceptable levels (peak-to-peak 
of few /iK). 

By focusing, as a working case, on the 100 GHz channels of Planck LFI, we have 
presented here extensive simulations of the straylight contamination starting from a wide 
set of simulated optical configurations, described in Paper I, in order to find the best 
compromise between resolution and GSC. 

By adopting some templates of Galactic foreground extrapolated from radio and IR 
surveys we found that it is possible to improve the angular resolution of about 5 — 7 % 
and to reach for example 10' — 12' of FWHM at 100 GHz by keeping the overall straylight 
contamination below the level of few /iK in terms of peak-to-peak and about 10 times 
smaller in terms of RMS, as necessary to avoid systematic errors comparable with the 
Planck sensitivity. 

We compared the level of straylight introduced by the different Galactic components 
for different beam regions and provided simple approximate relations giving the RMS 
and peak-to-peak levels of the GSC for the intermediate pattern and the far sidelobes as 
functions of the corresponding contributions to the integrated antenna pattern response, 
related to the edge taper. 

For the considered optical designs, the most important straylight contamination de- 
rives from the far sidelobes, where the Galactic signal overwhelms the other straylight 
contributions (nevertheless, in the intermediate pattern, the straylight contamination 
from CMB fluctuations is found to be not negligible compared to the GSC). 

We demonstrated that including the third order optical interactions changes only 
of few % the results of straylight analyses. As discussed in Paper I, this is extremely 
important for optical design optimization studies, probing that accurate enough optical 
simulations can be carried out by saving about 75% of the computational time without 
a relevant loss of accuracy. 

Keeping at very low level the GSC in the temperature data directly assures a reliable 
suppression of the GSC in the polarization data, at least in terms of RMS. The CMB 
polarization angular power spectrum (ET, E, and B modes) recovery will be then not 
significantly affected by straylight. 

The results at 100 GHz have been compared with those at 30 GHz, where the GSC is 
more critical, showing a peak-to-peak value at a level of ~ 5 — 7fiK. In comparison with 
previous analyses (see e.g. Burigana et al. 2001), this represents a relevant improvement 
in straylight rejection related to the optimization of the overall Planck optical scheme 
(Dubruel et al. 2000). Clearly, assuming different Galactic templates implies differences 
in the computed straylight signals. On the other hand, even for the quite different input 
maps adopted in these tests, differences larger than ~ 1 — 2/iK are limited to quite 
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localized regions, close to the Galactic plane, where the CMB anisotropy is dominated 
by the very high Galactic signal in the main beam. This suggests that, even far from 
the "cosmological window", a subtraction of the GSC well down to ~ 1/xK level could 
be obtained by evaluating with few iterations the straylight contamination signal for the 
maps directly derived from Planck observations, provided that the antenna pattern 
response could be quite accurately modelled. 

Acknowledgements. Some of the results in this paper have been derived using the HEALPix 
(Gorski et al. 1999). 

Appendix A: Transformation rules between telescope frame and beam frame 

Let s be the unit vector, choosen outward the Sun direction, of the spin axis direction 
and k that of the direction, z, of the telescope line of sight (LOS), pointing at an angle 
a ~ 85° from the direction of s. 

On the plane tangent to the celestial sphere in the direction of the LOS we choose 
two coordinates x and y, respectively defined by the unit vector i and j according to the 
convention that the unit vector i points always toward s and that x,y, z is a standard 
Cartesian frame, referred here as telescope frame. 

Let ibf,jbf,kbf be the unit vectors corresponding to the Cartesian axes Xbf,ybf,Zbf of 
the beam frame; k\ / defines the direction of the beam centre axis in the telescope frame. 
The beam frame is defined with respect to the telescope frame by three angles: 9b, (f>B, 4>b 
(8b and <j>B, two standard polar coordinates defining the direction of the beam centre 
axis, range respectively from 0°, for an on-axis beam, to some degrees, for LFI off-axis 
beams, and from 0° to 360°). 

Let i'bf, j' b f, k'bf {k'bf = kbf) be the unit vectors corresponding to the Cartesian 
axes x',y',z' of an intermediate frame, defined by the two angles 9b and 4>b, obtained 
by the telescope frame x, y, z when the unit vector of the axis z is rotated by an angle 

on the plane defined by the unit vector of the axis z and the unit vector kbf up to 
reach kbf. 

k'bf = hf = cos(^ B )sin(0 B )i + sin(^ B )sin(0 B )j + cos(9 B )k (A.l) 

i'bf = [cos(0 B ) 2 cos(6» B )+sin((/) B ) 2 ]2+[sin(0 B )cos((/) B )(cos(6' B )-l)]j-sin(6' B )cos((/) B )A:(A.2) 

j' bf = [sin((/) B )cos((/)B)(cos(6'B)-l)]«+[cos(6'B)sin(0s) 2 +cos(0B) 2 ]j-sm(6' B )sin(0B)fc .(A. 3) 

The beam frame is obtained from the intermediate frame through a further (anti- 
clockwise) rotation of an angle ipB (ranging from 0° to 360° 7 ) around kbf and is therefore 
cxplicitely given by: 

7 We note that, in other conventions, angles (j>' B and ip' B ranging from —180° to 180° are given, 
instead of cf) B and ip B . The angles 4>b and ipB here defined are equal to <j)' B and ip' B when they 
are positive and are given respectively by 360° + 4>'b an d 360° + ip' B for negative <f>' B and tjj' B . 
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i bf = [cos(i) B )i'bf,x + sm(ip B ) j'bf,x$ + [ c °s(ipB)i'bf,y + sm(tl> B )j' bfy ]j 
+[cos(ip B )i' b f,z + sm(ip B )j' bfz ]z 



(A.4) 



j bf = [sm(ip B )i' bfiX + cos(V>B)j' 6/; J« + [~sm(iP B )i'bf,y + cos(tP B )j' bfy ]j 

+[-sm(il> B )i'bf,z + cos(ip B ) j' bftZ ]z , (A.5) 

where the bottom index x (y, z) indicates the component of intermediate frame unit 
vector along the axis x (y, z) of the telescope frame, as defined by eqs. (A1-A3). 

Appendix B: Simulations at 30 GHz 

The feed horn at 30 GHz considered in these simulations is specified by its Spherical 
Wave Expansion (SWE) provided by Alcatel Space Industries, since the sub reflector 
is in the near field of the corrugated horn and near field effects cannot be neglected. 
The feed horn directivity is about 22 dBi, the ET is 30 dB at 22°, and the main beam 
has a FWHM resolution of 33.73'. The beam position and orientation is identified by 
(6 B , 4> B , i> B ) = (4.3466°, 153.6074°, 337.5°). 

We carry out the simulation of Planck observation as described in Sect. 2 by assum- 
ing both the Galactic templates presented in Sect. 2.2 and the WMAP map at 33 GHz. 
The results are summarized in Table 7 while Fig. B.l reports the TOD corresponding 
to the overall straylight signal (from the far sidelobes plus intermediate pattern) for the 
sum of the three Galactic components described in Sect. 2.2 and for the WMAP map. In 
spite of the differences in the foreground templates (the dominant signal deriving from 
Galactic diffuse free-free emission by using the templates of Sect. 2.2 and, according to 
Bennett et al. 2003b, from Galactic diffuse synchrotron emission by using the WMAP 
map) the peak-to-peak values is at a level of ~ 5 — 7/iK for the contributions from both 
the far sidelobes and the intermediate pattern (RMS ~ 1/^tK, mainly due to the signal in 
the far sidelobes). Right panel of Fig. B.l reports shows the difference between the TOD 
obtained by using these two different templates: only for ~ 0.01 % (0.09, 0.3, 10 %) of 
the samples of the TOD the difference is larger than 4/iK (3, 2, 1/iK). The figure shows 
also that differences larger than ~ 1 — 2/iK are localized quite close to the Galactic plane, 
where the CMB anisotropy is dominated by the very high Galactic signal in the main 
beam. This implies that a subtraction of the GSC well down to ~ 1/iK level does not 
require a particularly accurate description of the microwave sky emission nor particu- 
larly sophisticated computations, even at frequencies where the GSC is relevant, but it 
is mainly related to a good knowledge of the antenna pattern response. 

Clearly, the unsubtracted GSC is relevant at frequencies far from the "cosmological 
window". On the other hand, this results represent a significant optical improvement 
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Fig. B.l. The same as in Fig. 4, but for the simulation at 30 GHz and the overall 
straylight signal. For a better comparison the adopted temperature range is the same 
in the two panels, although the mimimum and maximum values are just different, as 
evident from the different peak-to-peak values reported in Table B.l. The right panel 
reports (antenna temperature in linear, not logarithm, scale, in this case) the difference 
between the middle panel and the left panel (see also the text). 

compared to the analyses of GSC at 30 GHz by Burigana et al. (2001), based on the 
optical simulations by De Maagt et al. (1998), predicting a similar level of GSC from the 
far sidelobes but a contamination from the intermediate pattern significantly worst (peak- 
to-peak of about 15/xK). This reduction of GSC from the intermediate pattern is due to 
the corresponding contribution to the integrated response from the antenna pattern at 
few degrees from the beam centre direction, significantly reduced in the actual optimized 
optical design 8 . 
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